Surface evolver simulations of drops on microposts 
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An important feature in the design of superhydrophobic surfaces is their robustness against collapse from the 
Cassie-Baxter configuration to the Wenzel state. Upon such a transition a surface loses its properties of low 
adhesion and friction. We describe how to adapt the Surface Evolver algorithm to predict the parameters and 
mechanism of the collapse transition on posts of arbitrary shape. In particular, contributions to the free energy 
evaluated over the solid-liquid surface are reduced to line integrals to give good convergence. The algorithm is 
validated for straight, vertical and inclined, posts. Numerical results for curved posts with a horizontal section 
at their ends show that these are more efficient in stabilising the Cassie state than straight posts, and identify 
whether the interface first depins from the post sides or the post tips. 

I. INTRODUTION 

A small liquid drop placed on the flat surface of a partially wetting material forms a spherical cap. A measure of the wettability 
of the surface is the contact angle 9y, defined as the angle between a tangent to the drop and the substrate, and related to the 
surface tensions between the three coexisting phases by Young's equation fH 

7sv — 7SL 

COS#Y = (1) 

7 

where 7sv, 7sl and 7 are the solid-vapour, solid-liquid and liquid-vapour surface tensions. Materials with contact angles 
9y > 90° or 6y < 90° with respect to water are referred to as hydrophobic or hydrophilic respectively. 

If a surface is covered by micron-scale bumps or posts, which can be stiff or flexible, its wetting properties are altered. In 
general, when a hydrophobic surface is roughened the contact angle is increased, leading to superhydrophobic behaviour J2)- 
Conversely, the contact angle of a hydrophilic material decreases as the surface becomes rough on a micron length scale Q. 
Micropatterned structures are relatively easy to fabricate using modern microlithography techniques, and superhydrophobic 
structures in particular are finding increasing applications as water repellent surfaces, energy efficient dehumidifiers or cooling 
devices[4|. A surprising number of plants and insects use superhydrophobic adaptations: examples include the water-repellent 
leaves of the lotus and nasturtium [5 1 and hairs on the legs of the water strider J6|. 

There are at least two possible states for a drop resting on a patterned surface. In the first of these, the collapsed or Wenzel 
state [7 1, shown in Fig.[TJa), the fluid fills the interstices between the posts. The effective contact angle in the collapsed state is 

cos 8w — r cos 6y (2) 

where the roughening of the substrate increases its surface area by a factor r. The Wenzel statel dominates the phase diagram 
for hydrophilic or weakly hydrophobic contact angles. 

The second common state is the suspended or Cassie-Baxter state [8 |; the drop lies on top of the posts with pockets of air 
beneath it, as shown in Fig.[TJb). The effective contact angle in the suspended state is 

cos 9c — 4> COS 9y — (1 — (j)) (3) 

where the projection of the liquid-solid contact area onto the base surface takes up an area fraction <f>. As expected the suspended 
state tends to be thermodynamically stable, or metastable, for larger By and more closely spaced posts, although if the posts 
overhang, it can be metastable even for hydrophilic contact angles ll9i fTTl . In the suspended state, a very high contact angle 
may be produced even for modestly hydrophobic materials and there is low resistance to droplet motion, an important property 
in many applications. Hence ensuring the stability and robustness of the suspended state is an important consideration in 
the application of both fabricated and natural rough surfaces, and it is interesting to calculate how this is affected by system 
parameters. 

For large drops, so that the interface under the drop is flat, and simple post geometries, the dependence of the collapse 
transition from the Cassie-Baxter to the Wenzel state on post dimensions and spacing and contact angle can be calculated 
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FIG. 1: (a) Collapsed and (b) suspended states on a superhydrophobic surface. 9w and 0c are the Wenzel and Cassie-Baxter angles respectively, 
and the "fa denote the surface tensions between adjacent phases. 




FIG. 2: Contact lines of various topologies (left), and the contours of their corresponding line integrals in (s, </>) (right), where the wetted area 
is enclosed by (a) a single contact line, not including the tip of the post, (b) a single contact line including the tip of the post, and (c) two 
separate contact lines wrapped around the post. 

analytically. However for more complicated surface geometries this quickly becomes impossible and numerical methods are 
required to understand how the instability to collapse occurs. Here we show how Surface Evolver, a public domain software 
package written by Brakke ifTZl \l3l which is widely used to predict equilibrium interface configurations, can be extended to 
identify the collapse transition on surfaces patterned by arrays of posts of arbitrary shape. 

Surface Evolver calculates the shape of a surface which minimises a given functional, in our case the free energy. The surface 
is parametrised in terms of triangular facets which are sequentially refined to give increasingly good approximations to the 
optimum surface configuration. Restrictions on the position of the surface are provided by boundaries, here the positions of 
contact lines between the surface and the posts. Hence to calculate the interface shape and free energy for a drop on an array of 
posts we need to define a suitable parametrisation for the position of the posts and to write the free energy contributions in terms 
of that parametrisation. 

We first describe how to parametrise post shapes, and then how to calculate and minimise the free energy of an interface 
in contact with the posts. Results are given for straight, inclined and bent posts; the first two of these can be compared with 
exact results to test the algorithm, the third can only be treated numerically. A discussion suggests possible applications of the 
approach. 



II. MODEL 



We consider thin cylindrical posts, with bases located at z = 0, arranged in a rectangular array with lattice vectors D x e x and 
D y e y . The posts are modelled by constructing a central spine of length L, described by the curve fq(s), where s is an arclength 
parameter, such that |r| = 1. The surface of the post is defined to lie at a specified perpendicular radius a(s) away from the 
spine. Hence it can be parameterised by s and an axial angle <f> as 

r(s, 4>) = r (s) + a(s) (— n(s) cos <j> + b(s) sin <f>) . (4) 

Here n is the unit normal to the spine in the direction of its curvature, such that 

r = kii , (5) 

where k is the scalar curvature, and 



b = r x n 



(6) 
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is the unit binormal (r'o, n, b form the Frenet-Serret frame). We assume posts to have a constant radius, except for a rounded 
end-section, included to avoid sharp edges for computational convenience, such that 

/ s J «o , < s < L - a , 

[ V a o ~ (L-a - s) z , L — a < s < L . 

We model the spine curve tq(s) of the post using a piecewise function described by a discrete data set. One choice is to model 
the spine as a series of straight line segments. However, this does not work well in conjunction with an interface, as it leads to 
sharp corners between segments which can pin the interface. A better method is to discretise the post using segments which are 
circular arcs, which for simplicity we take to all have unit length. Under this discretisation scheme, the data needed to describe 
apost of length N are the tangent vectors {r(p)} , p = 0, 1...N at each of the joins between segments. 

The curvature of each segment is 

k(s) = Kp, p +i = arccos [r(p).r(p + 1)] , (8) 

with p = [s\ (the largest integer smaller than s). The tangent vector within a segment is given by 

= r(p) sin [n(p + 1 - s)} + r(p + 1) sin [k(s - p)] ^ 
sin k 

which may be integrated to find the position of the spine, 

r (s) = r(p) + ^(P) cos i K (P + 1 ~ s )} + Kg + !) C 1 ~ cos l K ( s ~ P))) (10) 



KS1DK 



For a given segment, the binormal b is constant and is given by 



■ / \ i r(p)xr(p+l) 

b(s) = bp. p+1 = ; , (11) 

SUIK 

and the normal n is thus 

n(s) = b(s)xr(s) 

— r(p) cos [k(p + 1 — s)] + v{p + 1) cos [k(s — p)] (12) 
sin k 

We consider the array to be in contact with an interface spanning the x and y directions. For simplicity, we ignore edge effects 
and assume periodic boundaries so that the interfacial profile is the same for each post. The interface connects with the post at a 
contact line which, through Eqn. (ffl, can be described by a relation between s and <f>. 

IE. FREE ENERGY 

Next we write down the three contributions to the free energy of the system. The first is the interfacial energy between the 
water and air, 

7i = 7// dS. (13) 

The liquid-gas interface >4lg can easily be parametrised using facets within the Surface Evolver framework. An interfacial 
energy proportional to area is the default integral in the package, and the surface tension is chosen to be 1. Using the value 
7physicai = 0.072Nm _1 for water, and taking one unit of length in our simulation to correspond to 1 fim, the pressure in the 
corresponding physical system is given through the Laplace relation as 

A Preal = ^cai L ;"°° A Psim « 72kPaA Psim . (14) 

Tsimulation -^physical 

Since the micron scale is much smaller than the capillary length, any variation of Laplace pressure due to hydrostatics may be 
ignored. A suitable initial state is a horizontal interface, spanning the periodic boundaries in the x and y directions, and attached 
to the uppermost part of the hair. 
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The second contribution to the free energy is the wetting energy associated with the solid surfaces, here the surface of a post, 



wetting 



ff 1SL (r)dS+ [[ 1SG (r)dS 

JJAsl J J Abo (15) 

iff - cos 9 Y (r)dS+ ff ls G (r)dS , 

■jJAsi. JJAfn.UAfin 



where Als an d Alg are the areas of the substrate in contact with liquid and gas respectively. Young's law ([T]i has been used 
in the second step to write the substrate tensions in terms of the equilibrium contact angle 6y- Assuming that the total substrate 
area remains unchanged, the second term in the second line is a constant contribution to the free energy which may be ignored. 
Using Eqn. (HJ, the expression for a cylindrical element of area on a post surface can be written 

dS(s,cf)) = a {— dro(s) + (1 + an cos <j)) {— cos (f>n(s) + sin0b(s))} dsd(f> . (16) 

Thus, the free energy of the liquid wetting the post follows from Eqn. ( p"5j ) as 

.Fyy = — 7cos#y ff a\] a 2 + (1 + an cos <f>) 2 dsd<t> . (17) 



a is zero except at the tip of the post. If we restrict k to be small at the tip, then the expression ( 17 i simplifies to 



J" w w — 7cos#y // 1 + a 2 (1 + an cos </>) dsd<p , (18) 



and, using Eqn. (|7]), gives 



F w w — 7cos#y // &o (1 + a/«cos0) dsd(j) , (19) 



A 



which conveniently holds for both regions < s < L — a and L — a < s < L. 

Als is fixed rather than free and, as a result, it is not suitable to model this contact surface directly by facets, especially if the 
post is curved. This is because, for a hydrophobic substrate, there is a tendency for individual facets to grow or shrink in order 
to shortcut the substrate curvature, while for a hydrophilic substrate the facets may overlap one another or spread beyond the 
contact line. Thus the post is distorted away from its prescribed shape. 

Therefore, the integration over <j> is performed analytically, to give a line integral, 



Tw = — 7 cos 8y <p ao {4> + an sin <fi) ds . (20) 

This is evaluated numerically using the parameterised 'boundary' capability of Surface Evolver, with vertices specified in the 
parameter space of (s, eft) and then mapped to positions in Cartesian space by Eqn. ( ft). 

The boundary dA^s depends on the topology of the contact line, as shown in Fig.|2J In the case where the contact line loops 
around the trunk of the post (case (b) in the figure), there is a branch cut where cf> jumps from ir to — it and dA^s must continue 
to the end of the post at s = L in order to include the wetted area beyond the end of the contact line. It is also possible that 
states with two contact lines, when the post is fully submerged in some parts, may arise. These may be accounted for by 'gluing' 
adjacent edges together along the line <fi = tt, in which case the contour integral is taken along the path shown in Fig.[2jc). 

The third and final contribution to the free energy is the Laplace pressure, 

^Laplace = -Ap(VL-Vb) , (21) 

where Vl is the volume of the liquid phase, and Ap is the difference in pressure between the liquid and gas phases, and Vq 
is a constant. The Laplace pressure Ap is either variable, in which case the volume is constrained as Vl = Vq (J- is then the 
Helmholtz free energy), or fixed, in which case Vl varies (Gibbs free energy). In either case Ap induces a non-zero mean 
curvature of the surface. 

The Laplace energy may be rewritten as a surface integral over the boundary of the fluid, 

^Laplace = Ap ze z .dS . (22) 

J JAlgIJAls 
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FIG. 3: Schematic illustration of the hair geometries considered: (a) Straight hairs perpendicular to the base substrate, (b) Straight hairs 
inclined at an angle Q. (c) curved hairs with a section parallel to the base substrate, described by Eqn. \27\ . The diagrams indicate the position 
of the base of the drop in the Cassie state. 



Within Surface Evolver, the portion of ( 22 1 taken over Alg is most easily evaluated as a vector integral over the surface, but for 



the same technical reasons as discussed for the wetting integral (15 1, the contribution from the solid surfaces is most efficiently 
evaluated as a line integral. Therefore we write 



Jx = Ap I / ze z .dS 
'a lg 



+ / / (zq + a [—n z cos0 + b z sin0]) (— daz'o + a (1 + ancoscf)) [— cos0n z + sin</>6 z ]) d(f)ds 

Mls 



and perform the integral over s explicitly to give 

= Ap / / zdxdy + Ap d> I clzq [~n z sin <f> — b z cos (j> — azoq 



,2rl 



4> (n z + b 2 — kUzZq) + dz'o (n z sin cf> + b z 



+ j (2n 2 — kzq) cos 2(f> + | (n 2 z — b 2 z — kti z Zq) sin 2<fr] 

+ a 3 K [\b z n z (cos0+ \ cos 30) + \ (3n 2 z + bf) shx(f)+ ^ (n 2 z - bf) sm3(/)]^ds , (23) 
which is again integrated over the boundaries shown in Fig. [2] 



IV. RESULTS 



We begin by considering the threshold for collapse from the Cassie-Baxter to the Wenzel state for straight vertical posts, as 
depicted in Fig. |3ja). The shape of the interface when a Laplace pressure is applied will be of a nontrivial form, but for this 
particular geometry, it need not be known to calculate the critical pressure for collapse Ap max . This is because the system has a 
non-varying cross-section, and as such, an advance — Sz of the contact line will not change the shape of either the contact line or 
the interface. Therefore, the free energy change is ifTTl 

ST = {^2ira cos 6 Y + Ap (D x D y - na 2 ) } Sz , (24) 

and hence the interface will advance down the posts, leading to collapse, at a Laplace pressure 

— 27TO7 COS 9y 



Ap„ 



(25) 



Similarly, when the posts are straight but inclined at an angle £1 to the z direction (Fig [3jb)), a generalisation of Eqn. (25 i is 
easily derived by noting that the effect of inclination is to shorten the perpendicular distance between posts IfTTl 



A Pn 



—2narf cos 9y 

D x D y COS il — 7TC 



(26) 
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FIG. 4: Variation of depinning pressure with D x = D y = T + 2a for straight hairs, with theory (lines) and simulation results (markers) 
compared. Values are shown for the parmeters Oy = 105°, fl — 0° (blue solid line, * markers), 6y = 105°, Q = 45° (blue dashed line, + 
markers), 9y = 120°, Q = 0° (red solid line, x markers), and 9y = 120°, Q = 45° (red dashed line, o markers), 

These results provide a useful validation of the algorithm. In Fig.|4]we plot the variation of the Laplace pressure for collapse 
with post spacing for a square lattice of posts, with a = 1 and spacing S between posts such that D x = D y = T + 2a, for 
two representative contact angles. Each result is obtained by starting with the Laplace pressure at zero and increasing it in small 
increments; between each increase the surface configuration is iterated until equilibrium is reached. The final Laplace pressure 
at which the interface collapses, which corresponds to the limit of metastability of the suspended state, is then recorded. 

Fig.|4]shows that the Surface Evolver algorithm is able to reproduce the exact results well. As expected collapse is easier for 
less hydrophobic posts. Inclining the posts tends to stabilise the suspended state because the effective gap between the posts is 
reduced. 

We also present results for posts which are curved. The lower section bends in a circular arc, until it is pointing horizontally 
at an arclength sj, at which point there is a straight horizontal section, as shown in Fig.[3jc). This profile resembles the hairs 
found on the plastrons of certain aquatic arthropods lfT4l - [l6l , and posts of this shape can be fabricated artificially 1 18 -20 1. The 
corresponding spine function is 




sin 



2.sj 



< s < sj , 
sj < s < L . 



(27) 



We choose a = 1, L = 32 and sj — 16. In the y direction, we define the spacing between the posts, T, by D y = 2a + T. 
In the x direction, there is a greater degree of ambiguity in the definition of the spacing. We choose D x = T + L x — Sj + 2a. 
Fig.|6]shows results for the threshold pressure as a function of spacing for three values of Oy- The horizontal section of the posts 
tends to pin the interface (demonstrated in Fig.[5ja)), providing a means of stabilising the suspended state, even when the posts 
are not hydrophobic. 

The simulations also allow us to identify the path by which the interface moves down the posts, thus identifying the most 
vulnerable feature of a given post design. For smaller T and more hydrophobic dy, collapse occurs when the interface depins 
laterally, around the barrel of the hair, as shown in Fig. |3Jb). The threshold pressure for lateral depinning was calculated by 
Chrisp and Thorpe [ 17l to be 



Dl - (2asin6» Y ) - 2a cos ^ 



2 sin 9y 



cot 9 Y < _g_ 

COt9 Y > 2DZ 



(28) 



This result is compared to the simulations in Fig. [6] Eqn. (28 1 assumes a system which is uniform in the x direction (essentially a 
row of infinitely long hairs), which leads to the prediction being slightly higher than the simulation results. Fig.[5Jb) shows that, 



FIG. 5: Graphical output of the Surface Evolver simulation for curved hairs. The spines of the hairs are represented as red lines, and the 
interface is shown in blue: (a) a typical state where the interface is pinned to the horizontal sections of the hairs, (b) lateral depinning of the 
interface from the sides of the hairs, (c) tip depinning from the ends of the hairs. 




for finite hairs, lateral depinning occurs around the 'crook' of the hair at s = sj, not uniformly along the section sj < s < L as 
assumed in the derivation of Eqn. ( |28] l. 

When the hairs are not hydrophobic or when the spacing is sufficiently large, the interface depins from the tip, rather than the 
sides, of the hair, as depicted in Fig.|5jc). This occurs at significantly lower threshold pressure than Eqn. (28 1. For the cases 
Oy = 105° and 9y = 120°, a sharp drop is observed in the graph Fig.|6]as T is increased, marking the switching from lateral to 
tip depinning. 
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V. DISCUSSION 

An important feature in the design of superhydrophobic surfaces is their robustness against collapse from the suspended, 
Cassie-Baxter configuration to the collapsed, Wenzel state. Upon such a transition a surface loses its properties of low adhesion 
and friction. Here we have described how to extend the Surface Evolver algorithm to predict the parameters and mechanism of 
the collapse transition on posts of arbitrary shape. In particular, contributions to the free energy evaluated over the solid-liquid 
surface need to be reduced to line integrals to give good convergence. 

We have tested the algorithm against the exact results available for straight, vertical and inclined, posts. Numerical results for 
curved posts with a horizontal section at their ends show that these are more efficient in stabilising the Cassie state, and identify 
whether the interface first depins from the post sides or the post tips. 

Several insects have superhydrophobic adaptations to allow them to live at length scales where surface tension dominates. For 
example the legs of the water strider are covered with hairs to help it to float on the water surface and to leave the water easily as it 
jumps. The water spider has a body covered with stiff hairs to trap a bubble of air, the plasteron, which it uses forrespiration[21 1. 
Moreover, devices have been constructed which exploit meniscus-dominated buoyancy. These include microboats made from 
wire mesh[22 23], and a robot that mimics a water strider ll24l . The approach presented here will help to investigate the relation 
of the detailed shape of the hairs to biological function and hence in the design of such biomimetic devices. Indeed numerical 
modelling will be important in designing optimal micropatterned surfaces for a wide range of applications such as drag reduction, 
water harvesting and water repellency. 
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